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ABSTRACT 

Aims. The aim of this paper is to study the evolution of the broadband spectrum of one of the brightest and nearest quasars 3C 273. 
Methods. We analyze the data obtained during quasi-simultaneous INTEGRAL and XMM-Newton monitoring of the blazar 3C 273 
in 2003-2005 in the UV, X-ray and soft y-ray bands and study the results in the context of the long-term evolution of the source. 
Results. The 0.2-100 keV spectrum of the source is well fitted by a combination of a soft cut-off power law and a hard power law. 
No improvement of the fit is achieved if one replaces the soft cut-off power law by either a blackbody, or a disk reflection model. 
During the observation period the source has reached the historically softest state in the hard X-ray domain with a photon index 
r = 1.82 ± 0.01. Comparing our data with available archived X-ray data from previous years, we find a secular evolution of the source 
toward softer X-ray emission (the photon index has increased by Ar =s 0.3 - 0.4 over the last thirty years). We argue that existing 
theoretical models have to be significantly modified to account for the observed spectral evolution of the source. 
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1. Introduction 

3C 273 is a radio loud quasar, with a jet showing superlumi- 
nal motion, discovered at the very beginning of quasar research. 
Being one of the brightest and nearest (z=0. 158) qua sars, 3C 273 
was in tensively studied at different wavelengths (see lCourvoisierl 
d 19981) for a review). 

In the X-ray band 3C 273 is thought to be quite a unique 
source in the sense that in its X-ray spectrum the "blazar-like" 
emission from the pc-scale jet mixes with the "nuclear" compo- 
nent, which reveals itself through the presence of a large blue- 
bump and a so-called "soft excess" below ~ 1 keV. Both are 
typical for emission from the nuclei of Seyfert galaxies. During 
the observation period discussed in this paper the radio to mil- 
limeter component, usually associated with a jet, has been in its 
lowest state. Thus, our data provide an opportunity to study the 
"nuclear" component in more details and to test theoretical mod- 
els for the emission of 3C 273. 

Surprisingly, in spite of the low level of the "jet-like" 
contribution, our data do not reveal the presence of typical 
"Seyfert-like" features in the spectrum. In fact, the 0.2-100 keV 
spectrum is well fitted by a simple model composed of two 
power laws or a combination of a hard power law with a soft 
cut-off power law. We find that more sophisticated models, 
which include a black-body for description of the soft excess, 
or a disk reflection component give comparable or worse 
fits to the data, than the above simple model. Moreover, the 
"disentanglement" of Seyfert and jet components of the X-ray 



spectrum, as proposed by iGrandi & Palumbol d2004h . is not 
observed in our data set. 
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IWalter & Courvoisieti d 19921) have found an anticorrelation 
between the X-ray spectral slope and the ratio of X-ray to ul- 
traviolet (UV) luminosity of 3C 273 using quasi-simultaneous 
data from EXOSAT, GINGA and IUE obtained in the 80-s. They 
have interpreted this correlation as evidence that the X-ray emis- 
sion is due to thermal Comptonization of the UV flux. The pres- 
ence of optical monitors on-board of both XMM-Newton and 
INTEGRAL enables us to analyze correlations between the big 
blue bump vari ability and the variabi l ity in X-rays. Repeating 
the analysis of I Walter & Courvoisieti d 1992b with our data set 
we do not find the previously observed X-ray slope - UV/X-ray 
flux correlation. The disappearance of this correlation is even 
more surprising because of the appearance of a previously unob- 
served correlation between the UV and the X-ray flux from the 
system. Taking into account the fact that the typical variability 
time scales in the UV and the X-ray bands are different (months 
in UV and days in X-ray), true simultaneity of the observations 
could be important. Our XMM-Newton and INTEGRAL observa- 
tions are the first observations in which UV and X-ray data were 
taken at the same time. 

It is possible, however, that appearance/disappearance of 
the UV-X-ray correlation reflects a real evolution of the spec- 
tral state of the source during the last 30 years. Collecting the 
archival X-ray data for about 30 years of observations, we find 
that the source indeed evolves towards a softer X-ray state. We 
find that the source has reached its historically softest X-ray state 
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Table 1. INTEGRAL observations of 3C 273. 



Period 


T start 


Rev 


On Time ISGRI/SPI/JEM-X 








(ks) 


Exposure(ks) 


la 


2003-01-05 


28 


120.8 




lb 


2003-01-11 


30 


10.6 




lc 


2003-01-17 


32 


106.8 




1 








182.5/182.4/12.4 


2a 


2003-06-01 


77 


57.0 




2b 


2003-06-04 


78 


52.6 




2c 


2003-06-16 


82 


33.5 




2e 


2003-07-06 


89-90 


227.1 




2f 


2003-07-18 


93 


20.5 




2g 


2003-07-23 


94 


14.0 




2 








308.9/310.5/40.3 


3 


2004-01-01 148- 149 


26.4 


152.3/NA/13.8 


4 


2004-06-23 


207 


92.7 


67.1/NA/29.2 


5 a 


2005-05-28 320-321 


278.3 




5 b 


2005-07-09 


334 


102.5 




5 








269.1/367.2/68.7 



with photon index r = 1.82 + 0.01 in June 2003. This spectral 
slope is (by AF ^ 0.3 - 0.4) softer than the typical slope ob- 
served in the 70-s and 80-s. This is definitely larger than the 
typical short time scale fluctuations of the slope which are at a 
level of Ar^O.l. 

It is important to note that spectral evolution of the source 
should not be ignored in theoretical modelling of nuclear and 
jet emission from 3C 273. For example, the very hard spectral 
state o f the source in the late 70-s has led lLightman & Zdziarskil 
d!987l) to the conclusion that "compact source" type models 
dominated by pair production can not be applied to 3C 273 in 
spite of its high apparent compactness. However, this conclusion 
is not applicable for the present state of the source. Similarly, the 
most popular model for the hard X-ray spectrum of the source, 
inverse Compton emission from the pc-scale jet, is satisfactory 
for the present state of the source, but cannot explain the T ^ 1.4 
spectral index observed in 70-s. 

The paper is organized as follows: in Section 2 we present 
the sequences of observations and methods used for data reduc- 
tion and analysis. In section 3 we present the results obtained, 
and discuss them in Section 4. We then give a summary of our 
analysis in the last part of the paper. 



2. Observations and Data Analysis 

2.1. INTEGRAL observations 

Since the launch of INTEGRAL dWinkler et alj2003l) on October 
17, 2002, 3C 273 was regularly monitored. Details of the data 
used in the current analysis are given in Table[H The preliminar y 
analysis of the 2003 data was done bv lCourvoisier et all d2003al) . 
A detailed study of the multiwavelen gth spectrum of June 2004 
was presented bv lTiirler et alj d2006l) . 

During the INTEGRAL observations 3C 273 was all the time 
in a very low state, with a flux of several milliCrabs (see Section 
|3.2| >. To obtain better statistics for spectral analysis we have com- 
bined data taken within a few weeks, which resulted in five data 
sets. The combined exposure times for each data set and each 
instrument on-board of INTEGRAL are listed in Table Q] 

The data for all INTEGRAL instruments were analysed 
in a standard way with the latest version of the off-line 



Table 2. XMM-Newton observations of 3C 273. 



Period Date 


Revolution 


On Time 


PN Exposure Time 






(ks) 


(ks) 


1 2003-01-06 


563 


8.95 


5.97 


2 2003-07-07 


655 


58.56 


40.60 


3 2003-12-14 


735 


8.91 


5.77 


4 2004-06-30 


835 


62.91 


13.93 


5 2005-07-10 


1023 


28.11 


19.33 



scientific analysis (OSA ) package 5.1 distributed by ISDCQ 
dCourvoisier et al.ll2003bl) . 

Most of the time, observations were done in the 5x5 dither- 
ing mode. The staring mode was used only during revolutions 
30, 148, and 149. It is therefore not possible to use SPI data dur- 
ing the third period. There are also no SPI data during the fourth 
period, as the instrument was in the annealing mode during this 
observation. 

For the analysis of the INTEGRAL/ 'JEM-X observations we 
have selected science windows for which the source was less 
than 3.5° from the center of the field of view. During the first 
three periods JEM-X 2 was operating, while during the fourth 
and fifth period it was switched off, and data were taken with the 
JEM-X 1 telescope. 

To minimize the influence of the flares in the surrounding 
background of the charged particles we have not used pointings 
in which the detector countrate exceeds the mean of the given 
period by more than 5%. 

2.2. XMM-Newton observations 

The log of the XMM-Newton observations analyzed in this paper 
is presented in Table|2] Among all the observations of the source 
done with XMM-Newton we have selected those that were done 
while INTEGRAL was observing. 

The XMM-Newton Observation Data Files (ODFs) were ob- 
tained from the on-line Science Archival The data were pro- 
cessed and the event-lists filtered using xmmselect within the 
Science Analysis Software (sas), version 6.5 and the most up-to- 
date calibration files. All EPIC observations were performed in 
the Small Window Mode. In all observations the source was ob- 
served with the MOS1, MOS2 and the PN detectors. As the PN 
detector, in contrast to MOS 1 and MOS2, is not affected by pile- 
up in our observations we concentrate our analysis on the PN 
data only. Patterns 0-4 were used in the analysis. Background 
spectra were extracted from areas of the sky close to 3C 273 
showing no obvious sources. ARFGEN and RMFGEN were then 
run to produce the corresponding ARE and RMF files. 

Optical and ultraviolet fluxes were obtained with the optical 
monitor (OM) of XMM-Newton with omichain following the in- 
structions of the SAS Watchout pageQ and using the values pro- 
vided for an AGN spectral type. The resulting magnitudes are 
given in Table [3] No OM data for the selected filters are avail- 
able for Period 3. 



1 http://isdc.umge. ch/?Soft+download 

2 http://xmm.vilspa.esa.es/external/xmm_data_acc/xsa/index.shtml 

3 http://xmm.esac.esa.int/sas/documentation/watchout/uvflux.shtml 
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Fig. 1. top panel: 2003-2005 IBIS/ISGRI lightcurves of 3C 273in 20-60 keV (upper curve, black) and 60-100 keV (lower curve, 
red). The data are binned in 1-d intervals. For each time frame the date given at the bottom corresponds to the first data point, tick 
interval is 1 day. bottom panel: Hardness ratio F(6o-ioofeV)/F(20-60feV)' 



Table 3. Optical and UV magnitudes of 3C 273 during each ob- 
servational period. 





Period 1 


Period 2 


Period 4 


Period 5 


V 




12.53+0.01 


12.69+0.01 


12.57+0.01 


B 






12.95+0.01 


12.86+0.01 


U 




11.62+0.01 


11.76+0.01 


11.69+0.01 


UVW1 


11.32+0.01 


11.28+0.01 


11.44+0.01 


11.30+0.01 


UVM2 


11.34+0.01 




11.37+0.01 


11.19+0.01 


UVW2 




11.14+0.01 


11.36+0.01 


11.16+0.01 



3. Results 

3.1. Timing Analysis 

In the top panel of Figure Q] we show the 3C 273 IBIS/ISGRI 
lightcurves in the 20-60 keV (black points) and in the 60-100 
keV (red points) energy ranges. In this figure the data are binned 
in one day intervals. ISGRI data show that the hard X-ray emis- 
sion from the system is variable on at least a one day time scale. 
In order to investigate the nature of this variability we have pro- 
duced the "flux-flux" diagram for the 20-60 keV and 60-100 keV 
energy bands, shown in the bottom right panel of Fig. [8] Large 
error bars of the single-day flux measurements do not allow to 
see the general trend of the dependence between the 20-60 keV 
and 60-100 keV flux. However, rebinning the data by building 
the average of three adjacent (along the 20-60 keV flux axis) 
original data points in the flux-flux plot we find that the whole 
data set has a roughly linear behaviour (Spearman correlation 
coefficient is equal to R = 0.89, with a probability of random 
coincidence P nm d — 7 x 10~ 3 ). This indicates that the variability 




06.01.03 



07.07.03 



14.12.03 



30.06.04 



20.07.05 



Fig. 2. top panel: 2003-2005 XMM/PN lightcurves of 3C 273in 
0.2-1 keV (upper curve, black) and 3-10keV (lower curve, red). 
The data are binned in quarter of day intervals, bottom panel: 
Hardness ratio F (3 _i fev)/F(0.2-ifeV> 



of the source in the 20-100 keV band is mostly due to changes 
in the normalization, but not in the shape of the spectrum. 

Fig.|2]shows the XMM-Newton lightcurves in 0.2-1 keV (up- 
per curve, black) and 3-10 keV (lower curve, red) bands. One 
can see significant changes of the 0.2-1 keV flux from data set 
to data set, but the variability at the time scale of single obser- 
vation, ~ 10 ksec is negligible. The plot of the hardness ratio 
(bottom panel of Fig. [2]) shows that the 0.2-10 keV spectrum 
hardens toward the end of the observation period. 
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Energy, keV 



Fig. 3. Broadband spectra of 3C 273 for each period. Data are fit- 
ted with the simple phenomenological model described in Table 
[5] The new data are compared to the average flux (grey line) and 
the 1-cr variability range (grey area) based on historic data of 
3C273 dTurleret al.lll999i) . 



3.2. Spectral Analysis 

Our INTEGRAL and XMM-Newton observations provide a very 
broad band (5 decades in energy) set of quasi-simultaneous data 
which spans from the UV (~ 1 eV) up to the soft gamma- 
ray (~ 100 keV) band. The spectra of the source for all the 
five selected data sets are shown in Fig. [3] Taking into ac- 
count a wide range of theoretical models proposed for the X- 
ray spectrum of 3C 273 we adopt a "simple to complex" strat- 
egy in the spectral analysis and attempt first to find a "mini- 
mal model" which can provide a satisfactory representation of 
the data and then consider how changing and/or adding various 
theoretically motivated components improves or worsens the fit. 
The hydrogen column density was fixed to the galactic value 
N H = 1.79 x I0 20 cm- 2 ( iDickev & Lockmanll990l) . 



3.2.1. Phenomenological model 

Data in the 2-100 keV energy range are well fitted by a simple 
power law, see Table [4] In our analysis we have left the ISGRI, 
SPI and JEM-X intercalibration factors free to vary with respect 
to the PN camera. The significant variation of the JEM-X inter- 
calibration factor reflects a known change of the gain parameters 
of the instrument. It is possible to fit all the data with the same 
SPI (C = 1.41) and ISGRI (C, = 1.19) intercalibration factors, 
but we have preferred not to do this, as the XMM-Newton and 
INTEGRAL data are only quasi-simultaneous, and the different 
value of the flux measured by XMM-Newton and INTEGRAL can 
be due to the physical change of the source flux on a one month 
scale. The fit of the second data set can be improved by the ad- 
dition of a broad emission line as discussed below. 

The simplest satisfactory fit to the 0.2-100 keV data can be 
achieved with a model consisting of two power laws. A soft 
power law with a photon index F\ dominates in the energy band 



Table 4. Parameters of the single power-law model in 2-100keV 
energy range. The fit was done leaving free intercalibration fac- 
tors for ISGRI (C,), SPI (C s ) and JEM-X (C;) with respect to the 
PN camera. The parameter uncertainties correspond to 1 <x level. 

Period 1 Period 2 Period 3 Period 4 Period 5 



1.75+ 



Power law 

1 oi +0.01 i o a +0.01 i 7n+0.0 

1.81 _„„, l.»4_ m i-'O-o.o 
,8.79 + °-^ 10.10 +od4 

Intercalibration factors 



-o„0l 

.04 o rM+0.13 £ 79+0. 
0.16 iu - iw -0.07 °'" -0.12 -0.09 



Cj 
Q 

c s 
x 2 

(dof) 



i i 1 +0.09 
1-1 -0.09 

92+ 05 

-Q-Q5 

1 on+0.26 

l - Jy -0.25 



| Tc+0.05 
l. J -0.05 

i 99+O.O3 
1 • -0.03 
1 4Q+0.16 



i TQ+0.09 r\ an+0.11 
I.J -0.09 UJU „j| 

1.421°U °- 98 -> 



920 
(920) 



Statistics 
1673 950 
(1632) (1153) 



J -0.07 

1139 
(1205) 



1 66 +001 
-Q-Pi 
8 54 +m 

-0.09 

67 +0 05 

-Q-Q5 

1 90+0.04 

L67 -0.14 

1419 
(1390) 



* In units of 10 11 erg cm 2 s 




Energy, keV 



Fig. 4. A 3C 273 spectrum in the 4-7 keV energy band as ob- 
served with XMM '/PN in July 2003 (second data set). Data are 
fitted with a single power law. A broad iron line is marginally 
visible in the residuals (see text). 



below 2 keV and describes the so-called "soft excess". A hard 
power law with photon index T2 dominates in the hard X-ray 
band. An F-test shows that a slightly better fit to the data can 
be achieved if one uses a cut-off power law instead of the single 
power law for the description of the soft emission. 

Surprisingly, we find that already this simple model provides 
a very good fit to the combined INTEGRAL and XMM-Newton 
data in all the observations except the second one. In the sec- 
ond data set the addition of a broad redshifted emission line at 
Ef e = 6.65 + 0.07 keV in the rest frame of the object makes the 
fit acceptable. This line can be identified with a fluorescent K a 
iron line. A broad iron line was seen in p revious observations of 
3C 27 3 by different instruments, see e.g. lYaqoob & Serlemitsosl 
d2000h . The large width and marginal detection of the line (see 
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Table 5. Parameters of the phenomenological model fit to the 
data. When the uncertainties are not quoted, the value of the pa- 
rameter was fixed. 

Period 1 Period 2 Period 3 Period 4 Period 5 



r, 

F* 

r 0.2- 



F* 



E Fe 
cr Fe 
1** 



Cut-off power law describing Soft excess 

9 Af*+ 0M 9 m+0.08 9 37+O.I2 9 a,r+0.12 9 qo+0.10 
-Q16 -QQ8 -0.13 """-Q.13 ^■•'■'-0.11 

0.65« 0.52+ ; 0.73^0.65^ 0.58! g 

t £f|+1.08 9 i 1+0.16 9 Q/l+0.59 i qq+0.49 9 ^£+0.24 
O.UO 076 _„ 26 A.JH_ QAl l.O^_ 035 Z..JU _ 03g 

Second power law describing hard tail 



2 7Q+0.02 

sje+g 

6.65 
1 98+ 055 

I - 7o -0.39 

42 50+ 3232 



1.67 



+0.02 



1.73+ 



o A9+0.22 7ri+0.10 
0.10 °-" -0.17 '"-0.12 



1 09+O.OI 
^-O.OJ 

10.03 

Gaussian line 

6.65!°$ 6.65 6.65 

31 +o:o7 39+ 022 50+ 092 

"•-' 1 -0.06 v - 3y -0.l2 "' J "-0.50 

6 1 +L3 5 29 +3 31 3 14+ 3 91 



C, 



1 1 1 +0.09 
ip -0.09 

0.87+| 
1.32 



9 

-0.24 



Intercalibration factors 

1 O Q +0.05 

1.16!™ 1-38! 
1.43 



1 Q 1 +0.09 ri 07+0.10 

i.ji_ fl0? u.j/_ 8i9 

0.10 u - v/ -0.07 

-0.15 ~~ 

Statistics 



1 63 +a01 
-Q.fll 

o co +0.08 

°- J -o.io 
6.65 

9 tc+2.97 
-1.3Q 

17 14+ 5 - 83 
1 /.it 1231 

66+ 005 
119+8:^ 

1.62+H 



1305 2174 1360 1594 1832 
(1308) (2022) (1519) (1634) (1778) 
Statistics without iron line taken into account 
1316 2204 1364 1596 1836 
(df) (1310) (2025) (1521) (1636) (1780) 



X 

(df) 



* in 10 11 erg cm 2 s 1 units. 

*** in 10~ 5 photons cirT 2 s~' units. 



Figure |H do not allow us to firmly distinguish ionised and non- 
ionised iron. In all other spectra the addition of an iron line does 
not significantly improve the fit. For completeness we show the 
degradation of the fit if no line is included. The complete list of 
the fit parameters for our phenomenological model is given in 
Table H 



3.2.2. Modelling of the soft excess 

The "soft excess" below 1 keV is typical of Seyfert galaxies. 
Several competing models for the physical origin of this compo- 
nent are available. In particular, the most simple "phenomeno- 
logical" models of the soft excess include a (multi-temperature) 
black body component. However, the origin of such blackbody 
component is not clear, because the inferred temperature is usu- 
ally too high to be accounted for standard accretion disks. More 
complicated models include e.g. reflection of external X-ray 
emission from the disk. Attempts to fit the soft excess in our 
XMM-Newton data with different models give the following re- 
sults. 

A black body model for the soft excess does not provide an 
acceptable fit to the data. Best fit parameters for this model are 
listed in Table [6]The disk reflection componen t calculated us- 
ing the kdblur and ref lio n models of XSPEC (Ro ss & Fabianl 
120051: ICrummv et al .1120061) provides a better description of the 
soft excess, than the blackbody model. Among the parameters of 
the re £1 ion model we fit only the spectral slope of the power 
law emission of the central source, fluxes of the direct and re- 
flected components, and the ionization parameter of the gas £ 
= AttF/iih, where F is the illuminating energy flux, and n# 
is the hydrogen number density in the illuminated layer; in the 
ref lion model £ is restricted to 1 < ^ < 10000 erg cm s _1 ). 
Parameters of the model fit are given in Table|7] In this Table we 



Table 6. Parameters of the model which includes a black body 
to describe the soft excess. 

Period 1 Period2 Period 3 Period 4 Period 5 



T(eV) 



78.9+JJ 

1 OT+0.05 
1,7 -0.06 



87.2+™ 
l.OS+^j 



Black body 
79.0+^ 
1.21 



1 2-10 

E Fe 
<r Fe 
n" 

Fe 

Cj 
Q 

c s 



1 CQ9+0.007 



6.65 
0.10 
4 56 +1 64 

-1.57 



1 18+0.10 

la8 -0.1Q 
J3+0.06 



0.03 

Power law 

I O77+0.001 1 7O7+0.005 
18 ''-0.001 lJ °'-Q.QQ6 

9.88 + ^! 7.75 + ™ 



78.0+ 2 ? 

Q4+0.02 



6.72 



-0.01 '•'-'-0.04 

Gaussian line 
+0 04 6 65 

0.06 U - UJ 

0.10 0.10 
4 89+ - 37 3 54+ 1 - 29 
Intercalibration factors 



i orc+0.005 

6 52 +t> " 2 

°- JZ '-0.03 



6.65 
0.10 

7 4Q+0.94 
^•^-0.95 



1 

1.72 



§1 

-0.31 



i -in+0.05 i 41+O.IO 
" J -0.05 1 -^ 1 -0.10 

1.34+!!" 3 1.80+0 1 :! 2 



1.65 



40 +0 12 

1 9 i +0.08 
^ -0.08 



82.0!^ 

-0.03 



1 71Q+0.004 
L - ' 1 ^-0.003 
o ^9+0.03 
-0.02 



6.65 
0.10 

9 QC+0.90 
i.JJ_ Q85 



70 +0 06 

' w -0.06 



-0.18 

Statistics 



1.98+ 



x 2 


1445 


2659 


1543 


1799 2161 


(dof) 


(1310) 


(2024) 


(1515) 


(1630) (1780) 




Statistics without iron line taken into account 


Xnl 


1453 


2710 


1551 


1807 2169 


(dof) 


(1311) 


(2026) 


(1516) 


(1631) (1781) 



In units of 10" 11 
* In units of 10 



erg cm ' s 

5 nl.nfn.. W\7~l 



photons keV cm s 



Table 7. Parameters of the model consisting of emission from 
the central source reflected from the relativistically-blurred pho- 
toionized disc and a jet 

Period 1 Period 2 Period 3 Period 4 Period 5 
Direct and reflected emission from the central source 



F* 

r 0.2-1 

1*1 

R(%) 



E Fe 
"1- 

Cj 

Cj 
c s 

x 2 

(dof) 

Xnl 

(dof) 



2977+1565 

-0 72 
9 oq+o:07 
z "° -0.08 
3.98 
1.00 



10000+° 

4.3 1! - 31 
2 56+^ 

z " JU -0.02 

20 ^ 02 +i79.5i 



3170!, 33 g 10000!° 
a sq+o.60 a ne+orra 

-0.59 w -0.39 

2 48+ 007 

z "^°-0.09 
10.69 



2 56+" 08 

-0.08 

2.61+ 201 



230+3.98 20 02 +1/9: " 2 61+ 2 " 1 1138 

^■^-Xm ■ iu - u ' i -2.39 u -1.56 il - JO -7.11 

Power law describing let emission 
1.53+™ 4 1.55+ 005 
6 



1-70!^ 
8.05!«| 



I 7C+O.OI 
1- ' J -Q.02 
q 99+0.25 
7 ' J -0.36 



6.65 

77+0.09 

1 16-8:2? 

I - IO -0.21 



-0.75 

Gaussian line 



-0.06 

5 76 +2 63 



6.65 



2.28'| 



3.06!» 6 62 

Intercalibration factors 
J9+0.02 j 28 +005 



6.65 

20 +0 W 

"•^"-0.20 



1. 

1.17! 

1.30! . 



1.12 



-0.09 



35 +0 05 
u -' u -o.io 



14 

Statistics 



2570 +23 » 1 



2 76 +() 04 
^•'"-o.oi 

2.69 +34 i 



■ <r 4 +o.03 

i - J ^-0.02 

8.02il 



6.65 
1 49+0.94 

i '^ z '-0.86 

0.66!^ 
1.08 + ° 
135!] 



1399 2239 1349 1516 1819 
(1410) (2002) (1507) (1542) (1758) 
Statistics without iron line taken into account 
1405 2262 1353 1516 1822 
(1412) (2004) (1509) (1543) (1759) 



+ in units of erg cm s~ 
* in units of 10" 11 
** in units of 10~ 4 photons keV" 1 



erg cm 2 s 1 . 



also report the measured ratio R of the reflected flux to the total 
flux expressed in percents. 

One should note, however, that the amount of reflection 
found in each data set, except for data set number 2, is at the 
level of several percent. Thus, a better fit (compared to black- 
body) to the data is explained not by the careful account of re- 
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flection features, but by a large contribution of the direct power 
law illuminating the disk, so that the disk reflection model is 
not significantly different from the previous phenomenological 
model. To summarize, we find that neither the black body nor 
the disk reflection provide better descriptions of the soft excess 
than the (cut-off) power law. 

Compton reflection is characterized by a "reflection hump" 
at about 30 keV. To look for the presence of a reflection compo- 
nent in the spectrum we have also tried to fit the spectra above 
2 keV with either a simple power law or a direct plus reflected 
power law (the re £1 ion model of XSPEC) and compare the 
quality of the fits. Although we find that the addition of the re- 
flected component at the level of about 20% improves the quality 
of the fits, the F-test shows that the probability that this improve- 
ment is coincidental is at the level of 10%, so that no definitive 
conclusion about the presence of a reflection hump in the spec- 
trum can be made. One should note that the largest uncertainty is 
due to the fact that the ISGRI instrument is sensitive only above 
~ 20 keV, while XMM-Newton data are available only below 
~ 10 keV and, in fact, there is a degeneracy between the normal- 
ization of the reflection hump and the ISGRI - XMM-Newton 
intercalibration factor. 

A combination of black-body and reflection components, 
suggested as a best-fit model for the Beppo-SAX data by 
IGrandi & Palumbol d2004l) does not change the above conclusion 
(see below), leaving the (cut-off) power law to be the best model 
for the soft excess. 



Table 8. Black body, power law, Seyfert-like component and 
Gaussian line fit to PN, JEM-X, SPI and ISGRI 3C 273 spectra. 
The fit was done leaving free intercalibration factors for ISGRI 
(Q), SPI (C s ) and JEM-X (Cj) with respect to the PN camera. 
The parameter uncertainties correspond to 1 <x level. 

Period 1 Period 2 Period 3 Period 4 Period 5 



f; i.oi!°- 



F » c 95+0.09 
2 J - z ~'-0.09 



p* t. 47+0.04 

n" 6 45 +L82 

n Fe -3.58 

Cj 1.09*» 

Ci 0.68! | 

c i or 18 

A' 1 • v - / * _Q j 9 



X 1340 
(dof) (1311) 

Statistics 

X 2 nl 1345 
(dof) (1312) 



Black body 



< 0.004 



53+ 03 
Power law 

f. qt+0.03 c Q^+0.07 

Seyfert component 

■2 Q£+0.01 T in+003 
^• s<:) -o.oi 2 - 1U -0.03 
Gaussian line 

3.03 + ?-£ 4.45 +274 



33+ 02 



4 76+ 04 

^•'"-0.05 



1 96 +0 02 

l y -0.02 



1.78 

Intercalibration factors 



+ 1.15 
1.78 



79+ 03 

u - ' '-0.02 



6 74+ 05 



i 7 1 +0.02 
1 -' 1 -0.02 



2.92 + 



1 QQ+0.10 

1-J -o.io 
1 16 +0W 



0.28!»;» 



65 +0 05 

1 05 +m 

1 -" J -0.02 
i 49+0.12 
1 -^ z '-0.12 



1 26 +om 

0.80!ol 
0.97!«ij 

Statistics 

2342 1383 1561 1913 

(2026) (1531) (1565) (1781) 
without iron line taken into account 

2350 1389 1562 1918 

(2027) (1532) (1566) (1782) 



* in units of 10 11 erg cm 2 s 

** in units of 10~ 5 photons keV~' cm -2 s~ 



3.2.3. Disentangling the Jet and Seyfert contributions 



IGrandi & Palumbol d2004l) (GP) have put forward a strong claim 
that Seyfert and "jet" components of the X-ray spectrum could 
be disentangled in the case of 3C 273 so that the observed vari- 
ations of the X-ray spectral shape are successfully described 
by variations in the normalizations of the template "jet" and 
"Seyfert" components. We have tested th is hypothesis with 
our da ta. For this we repeat the analysis of IGrandi & Palumbol 
d2004l) . First, we find the "template" models for Seyfert and jet 
components by fitting simultaneously all the 5 data sets with the 
same model which consists of a hard power law (jet) and black- 
body plus partially reflected cut-off power law (Seyfert compo- 
nent). The best fit is achieved with the following parameters: a 
black body temperature Tj,b = 75 eV, a Seyfert component power 
law slope r s = 2.25, a reflection component fraction RC = 0.77, 
a jet power law slope T y = 1.5, an inclination of the disk <p = 18°. 
We assumed solar abundances of elements and fixed the posi- 
tion and width of the iron line at the values found in the second 
data set: Ep e — 6.65, o- Fe = 0.3 keV. Then we froze the above 
parameters and left only the normalizations of the black body, 
reflection and jet power law components variable. With such a 
"restricted model" we fit each data set separately and find the 
normalizations of the three "template" components. The fluxes 
in each component found in this way are listed in Table[8] 



After this analysis, and following lGrandi & Palumbol d2004l) . 
we study the correlations between e.g. Seyfert component and 
iron line strength or between fluxes in Seyfert and black body 
component, which were found in the Beppo-SAX data. Fig. [5] 
(which is an analog of Fig. 1 of GP) shows that no such correla- 
tions are present in our data set. We discuss the implications of 
this result below. 



3.2.4. Multi band flux and spectral slope correlations 

Comptonization of big blue bump photons was proposed as a 
plausible mecha nism of production of the har d X-ray emission 
from 3C 273 by IWalter & Courvoisierl dl992l) . The motivation 
for this model was a linear correlation between the X-ray spec- 
tral slope above 2 keV and the logarithm of the ratio of the 2- 
lOkeV flux to the ultraviolet and soft excess count rates found in 
EXOSAT, GINGA and WE data. The presence of the UV camera 
on board of XMM-Newton and of the Optical Monitor Camera 
(OMC) on board of INTEGRAL enables us to test this correla- 
tion between the UV and X-ray fluxes and slopes. Fig. [6] shows 
the values of the power law index in 2-10keV band, Y2, plot- 
ted against the logarithm of the ratio of X-ray and UV fluxes 

(F2-l(Uev/vFi0cv)-U 

Figure [6] clearly shows that the correlations seen in 
EXOSAT/GINGA X-ray data and IUE UV data are not present 
in our data set. For completeness we have extended the anal- 
ysis for the whole set of publicly available XMM-Newton data 
on 3C 273 to produce this Figure (for details see Soldi et al., in 
preparation). 

In Fig. [6] two important facts are apparent. First, the scatter 
of the XMM-Newton data points along the x-axis is much smaller 
than that of the EXOSAT/GINGA and IUE points. This means 
that in the recent past the X-ray and UV fluxes have varied in a 
correlated way. Left top panel of Fig.[8]shows that this is indeed 
the case (Spearman correlation coefficient is equal to R — 0.93, 
with a probability of random coincidence P mn d — 3 x 10~ 3 ). This 
is a quite surprising result, because previous searches for such 
correlations gave negative results. However another correlation 
typical for Seyfert galaxi es between the hard X -ray spectral in- 
dex F2 and the X-ray flux dZdziarski et alj d2003t) . and references 



4 We have extrapolated the UV spectrum obtained with the XMM- 
Newton optical monitor up to 10 eV using the time-averaged UV spec- 
trum of 3C 273. 
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Fig. 5. The dependance of the 0.2-2 keV black-body emission (left), the 2-10 keV jet like flux (middle), and the iron line flux (right) 
on the 2-10 keV Seyfert-like flux. Data points are taken from Table [8] 
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Fig. 6. 2-10 keV spectral photon index of the phemenologi- 
cal model as a function of the logarithm of the ratio of the 2- 
lOkeV to ultraviolet fluxes. Old Ginga and EXOSAT data are 
shown w ith red triangles. The l i ne ind icates the best correlation 
found in lWalter & Courvoisierl d 1992b . New XMM-Newton data 
are shown with black circles. We have extrapolated the UV spec- 
trum obtained with the XMM optical monitor up to 10 eV using 
the time-averaged UV spectrum of 3C 273. 



GINGA and IUE data and "appearance" of X-ray - UV flux cor- 
relations not detected before. 

Having found a previously undetected correlation between 
UV and X-ray flux, we have searched for the presence of addi- 
tional correlations, such as between the 0.2-1 keV and 3-10 keV 
fluxes (top right panel of Fig.O, 0.2-1 keV and 20-60 keV fluxes 
(bottom left panel of Fig. [8]> and 20-60 keV and 60-100 keV 
fluxes (bottom right panel of Fig. [8}. One can see a clear correla- 
tion between 20-60 keV and 60-100 keV fluxes, indicating that 
both bands are dominated by one and the same physical compo- 
nent (e.g. powerlaw emission from the jet). On the other side, no 
or only marginal correlation in the two other cases shows that 
several physical mechanisms contribute to the emission in these 
energy bands. 



4. Long-term evolution of the X-ray spectrum of 
3C 273 

To test the evolution of the X-ray spectrum over the last 30 year 
period we have collected the archival X-ray data from different 
missions. The evolution of the hard (2-10 keV) s pectral index 
found from the archival dat£0 ( Tiirler et al.l 1 1999b is shown in 
Fig. [9] One clearly sees both significant short-term variations 
of the spectral index with a scatter of Ar ^0.1 and a long- 
term trend of softening of the spectrum by as much as AT ^ 0.4 
over the last 30 yr. It is interesting to note that during our 2003- 
2005 monitoring campaign the source has reached its histor- 
ically softest hard X-ray spectrum characterized by a photon 
index F = 1.82 ± 0.01. It is possible that still softer spectral 
states were found at the very beginning of X-ray observations 
(UHURU, ARIEL5/C and OS08/A data points), but the uncer- 
tanity of those measurements are too large to be conclusive. 



th erein), is no t obser ved (see Fig . O . A similar result was found 
bv lPageetal]d2004 . 

Second, the values of the hard X-ray spectral index in XMM- 
Newton observations are systematically higher than those of the 
EXOSAT/GINGA data. This indicates that the spectral state of 
the source has evolved over the 30 years of X-ray observations 
(in the next section we discuss this in details). This may be linked 
to the "disappearance" of the correlations seen in EXOSAT, 



5. Discussion 

Most of the theoretical models of 3C 273 suggest that the source 
spectrum contains two major contributions, one from a pc-scale 
jet and the other from a Seyfert-like nucleus. The very low mil- 
limeter flux measured during our observations suggested that 
the jet contribution has reached its historical minimum (see 

5 Available from 3C 273's Database: |http://isdc.unige.ch/3c273/] 
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Fig. 8. Flux-flux plots for different energy bands. On the right bottom panel one day averaged INTEGRAL data are shown. To reduce 
the scatter the data were rebinned by building the average of three adjacent original data points in x direction. 



iTiirler et"aT] d2006l) ). However, this did not result in the clear ap- 
pearance of the Seyfert-like contribution. This means either that 
the Seyfert-like contribution has decreased simultaneously with 
the jet contribution (which, because of the difference in the size 
scales and production sites can happen only by a chance coinci- 
dence), or that the observed emission just can not be decomposed 
into Seyfert and jet contributions. 

Our data show that the separation of Seyfert and jet contri- 
butions to the spectrum could not be done by fixing a "template" 
for each of them. Namely, although the typical ranges of vari- 
ations of soft and hard X-ray fluxes in our data set are simi- 
lar to the ones of Beppo-SAX data set, the "template" models 
for the Seyfert, jet and black body components found in our 
data set are significantly different. In particular, the additional 
black body component is not required by the data and thus the 
resulting black body flux shown in the left p anel of Fig.[5]is sys- 
temati cally lower than the values reported in lGrandi & Palumbol 
d2004l) . Our data show that the correlation between the black 
body and Seyfert components is absent at such low black body 
flux values. The same is true for the non-observation of the cor- 
relation between the Seyfert component flux and the strength of 
the iron line: since the addition of the iron line is not statistically 
required in any data sets except for the second one, no corre- 



lation is seen at such low iron line flux values. However, it is 
important to note that no direct compa rison between our Fig. [5] 
and Fig. 1 of iGrandi & Palumbol d2004l) can be done, because of 
the differences in the template models derived from Beppo-SAX 
and INTEGRAL -XMM-Newton data sets. 

In fact, the analysis of the secular evolution of the source 
presented above shows that both Seyfert and jet spectra evolve 
with time. This can be illustrated by the results presented in Fig. 
[6] where we plot the hard X-ray spectral slope as a function of 
the ratio of the X-ray power law flux to the UV flux. This ratio 
is apart from a constant factor the ratio of the Comptonized flux 
to the source photon flux in a model in which the X-ray power 
law is due to the Comptonization of the soft photons by a ther- 
mal electron population. The points given by the red crosses with 
large error bars were obtained by t he EXOSAT and IUE m i ssions 
in the 80-s and were analysed by Walter & Courvoisierl d 19921) 
who deduced that in the frame of the thermal Comptonization 
model, the electron temperature was about lMeV or some- 
what less. They also deduced that the optical depth of the 
Comptonizing medium was 0.1-0.2 and that it covered a cou- 
ple of percent of the soft photon source. 

The XMM-Newton data shown in Figure [6] by the black 
points with small error bars clearly do not match the correlation 
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Fig. 9. Evolution of the flux density at 1 keV, as extrapolated from the fit done above 2 keV, in lO^photons/keV/cm 2 /.? (top panel) 
and hard X-ray spectral index (bottom panel) over the 30 years of observ ations. UHURU, Ariel 5/C, OSO 8/ A, HEA O-1, HEAO- 
2/MPC , and BALL/KIT dat a points are taken from iMalaguti et al.l dl994l) : EXOSAT and Ginga data are taken from lTurner et al.l 
dl990l) : ASCA data are from lCappi et all (Il998l) : RXTE, Beppo-SAX, and part of XMM-Newton data are from Soldi et al. (in prepa- 
ration). 



found in earlier years. This is particularly true for the spectra ob- 
tained during periods 1 and 2, when the soft X-ray emission was 
much brighter than at the other epochs, while the ultraviolet flux 
was similar. This indicates a much harder blue bump for these 
epochs. More photons were therefore available to Compton cool 
the hot electrons than inferred from the ultraviolet flux shown 
in Figure 5. This therefore leads to the softening in the X-ray 
emission. Alternatively, one can also infer a modification of the 



Comptonization region when the source is very weak. The re- 
cent data also show that the hard X-ray flux is correlated with the 
soft flux (see Figure 8), an effect not observed at earlier epochs 
(ICourvoisier et al.lll990T) . 

The correlation between the optical/UV and the X-ray flux 
is characteristic for Seyfert galaxies. This relationsh ip is particu- 
larly c lear in the long-term behaviour of NGC 5548 dUttlev et al.l 
120031) . but was also seen in other Seyfert galaxies like for 
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Fig. 7. The dependence of the hard X-ray photon index (derived 
from fitting the data in 2-100 keV energy range with a single 
power law, see table|4ji on the 2-10keV X-ray flux. 



instance NGC 4051 dShemmer et al.l d2003l) . and references 
therein). The data of 3C 273 presented here do also show the 
appearance of such a trend (see Fig. [8j, apparently suggesting 
a Seyfert-like behavior of 3C 273 in 2003-2005. However, it is 
clear from the above discussion that in the case of 3C 273 con- 
clusions of this sort should be taken with caution: there is no 
clear way to disentangle the Seyfert and jet contributions. For 
example, another characteristic of Seyfert galaxies, a correla- 
tion between the hard X -ray spectral index F2 and the X-ray flux 
dZdziarski et al.l d2003l) . and references therein), is not observed 
(see Fig. |7]i. 

The secular evolution of the hard X-ray spectral index over 
the last 30 years, clearly seen in Fig. [9] shows that the jet con- 
tribution also evolves with time. Whether the hard X-ray flux is 
generated in the jet, or it is due to the nuclear related physic s 
(i.e. from the "compact source" dLightman & Zdziarskilll987D ) 
is still completely undecided. "Jet" and "compact source" mod- 
els differ in their approach to the so-called compactness problem. 
The essence of the compactness problem is that the naively es- 
timated optical depth of the source for y-rays with energies E y , 
such that pairs are produced in the interactions with soft pho- 
tons of energies e s > 2m 2 e c A /E y , is higher than 1. The "compact 
source" models self-consistently a ccount for the effect of the pai r 
production cascade in the source dLightman & Zdziarskilll987l) . 
The "jet" models get rid of the compactness problem by assum- 
ing relativistic beaming of the y-ray emission. 

The observed evolution of the hard X-ray spectral index from 
r =* 1.4 to F =a 1.8 over the last 30 years could help to distin- 
guish between the two possible models of the high-energy emis- 
sion. The point is that although the present-day values F > 1 .5 
are easily explained in both models, the very hard spectrum ob- 
served is the late 70-s is not. In both models the gamma-rays are 
produced via inverse Compton emission from high-energy elec- 
trons. The hard spectral index F ss 1.4 observed in 70-s and 80-s 
(see Fig. |9]l is too hard for the optically thin inverse Compton 



emission from an electron distribution cooled via synchrotron 
and inverse Compton energy loss which can produce only spec- 
tra with T > 1.5. 

This problem can, however, be relaxed if one assumes that in 
addition to electrons the emission region contains also protons. 
If the main source of the seed photons for the inverse Compton 
scattering is UV photons from the big blue bump, the hard X-ray 
inverse Compton emission originates from a population of elec- 
trons with moderate energies E e ~ 1-10 MeV. In this energy 
range cooling could be dominated by ionization, rather than in- 
verse Compton/synchrotron losses. The Coulomb loss leads to 
the hardening of the electron spectrum, compared to the injec- 
tion spectrum by AT„ ^ - 1 . This will result in a hardening of 
the inverse Com pton spectrum be low the so-called "Coulomb 
break" (see e.g. lAharonianl d2004l) ). A simple estimate shows 



that if the proton density in the gamma-ray emission region is 



10 10 cm 3 , the Coulomb break is situated in the hard X- 



ray band which can provide the explanation of the harder than 
r = 1.5 spectrum of the source in late 70-s. High matter density 
in the emission region obviously favors the "compact source" 
models compared to the jet models: the assumption about a large 
(10 cm 4 ) matter densities at the pc distances from the nucleus 
needed to explain the hard spectrum in the jet model framework 
looks somewhat extreme (for comparison, the density of dust 
and molecular gas in the Central Molecular Zone (CMZ) of the 
Galaxy is estimated to be 10 4 cm 3 , see e.g. lMorris & Serabvnl 
d 1996b '). The conclusion that the observed evolution of the hard 
X-ray spectral index favors the compact source models should 
also be taken with caution. The reason for this is that the pres- 
ence of a Seyfert-like contribution to the hard X-ray spectrum 
can affect the value of the spectral index. 

If the gamma-ray emission of 3C 273 originates in a compact 
source the synchrotron radiation from electrons which upscatter 
the blue bump photons to the GeV energy band could contribute 
also to the soft X-ray spectrum of the source. The magnetic field 
strength in the y-ray emission region is poorly constra ined by 
the observations (see, however. ICourvoisier et alj d!988l) for re- 
strictions on the magnetic field strength during flaring activity of 
the source in the late 80-s). Assuming that the energy density of 
magnetic field in the compact emission region is comparable to 
the energy density of radiation, 



U 



P h 



2xlO z 



10 46 erg/s 



10 16 cm 



R 



erg/cm 



(1) 



(we estimate the size of emission region to be R ~ 10 16 cm from 
the observed intra-day variability of the source) one obtains an 
estimate 



L 


1/2 


10 16 cm 


10 46 erg/s 




R 



G 



(2) 



The synchrotron radiation from 1-10 GeV electrons in such mag- 
netic field is emitted at the energies 



B 



10 2 G 



10 GeV 



keV 



(3) 



The n on-observation of a cut-off in the GeV band dCollmar et al.l 
l2000l) implies that the high-energy cut-off in the electron spec- 
trum is at least at energies E e > 1 GeV which, in turn, implies 
that the cut-off in the synchrotron spectrum should be at least in 
the soft X-ray band. Taking into account that the best fit to the 
"soft X-ray excess" in our data is given by a simple cut-off power 
law model, we find that it is possible that the observed soft X-ray 
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Fig. 10. Synchrotron-external-Compton fits to the data for the 
observation period 2 in the "compact source" type model. Solid 
line: the total spectrum, short dashed line: the inverse Compton 
spectrum without the account of the pair production. Dotted line: 
synchrotron spectrum. Long dashed line: the seed photon spec- 
trum for the inverse Compton emission. Similar to Fig. [3] the 
average flux (grey line) and the 1-cr variability range (grey area) 
is shown. 



excess in the spectrum of 3C 273 is the high-energy tail of the 
synchrotron radiation from the compact core of the source. As 
an example, we show in Fig.[l0]a synchrotron-inverse Compton 
fit to the broad band spectrum of the source for the observation 
period 2 (with the historically steepest spectrum in the hard X- 
ray band). The parameters used for the model fit are: a magnetic 
field B = 33 G, a size of emission region R = 10 16 cm, a spec- 
tral index and cut-off in the electron spectrum are T e = 2.6 and 
E cu t - 10 GeV, respectively. The gamma-ray spectrum has a 
break at £ = 10 MeV due to pair production. For comparison, 
short-dashed line in this figure shows the spectrum of inverse 
Compton emission without taking into account the pair produc- 
tion. 



lar evolution. The secular evolution of the source is evident from 
the study of the evolution of the X-ray spectral index during the 
last 30 years of observations (Fig. [9]). We argue that existing the- 
oretical models have to be significantly modified to account for 
the observed spectral evolution of the source. 



7. Acknowledgments 

Authors are grateful to the unknown referee for helpful 
comments. PL acknowledges the support of KBN grant 
1P03D01128. 



References 

Aharonian, F. A. 2004, Very high energy cosmic gamma radiation: a crucial win- 
dow on the extreme Universe (published by World Scientific Publishing) 
Cappi, M., Matsuoka, M., Otani, C, & Leighly, K. M. 1998, PASJ, 50, 213 
Collmar, W., Reimer, O., Bennett, K., et al. 2000, A&A, 354, 513 
Courvoisier, T. J.-L. 1998, A&A Rev., 9, 1 

Courvoisier, T. J.-L., Beckmann, V., Bourban, G., et al. 2003a, A&A, 411, L343 
Courvoisier, T. J. L., Robson, E. I., Blecha, A., et al. 1990, A&A, 234, 73 
Courvoisier, T. J.-L., Robson, E. I, Hughes, D. H., et al. 1988, Nature, 335, 330 
t> Courvoisier, T. J.-L., Walter, R., Beckmann, V., et al. 2003b, A&A, 41 1, L53 
Crummy, J., Fabian, A. C, Gallo, L., & Ross, R. R. 2006, MNRAS, 365, 1067 
Dickey, J. M. & Lockman, F. J. 1990, ARA&A, 28, 215 
Grandi, P. & Palumbo, G. G. C. 2004, Science, 306, 998 
Lightman, A. P. & Zdziarski, A. A. 1987, ApJ, 319, 643 
Malaguti, G, Bassani, L., & Caroli, E. 1994, ApJS, 94, 517 
Morris, M. & Serabyn, E. 1996, ARA&A, 34, 645 
Page, K. L., Turner, M. J. L., Done, C, et al. 2004, MNRAS, 349, 57 
Ross, R. R. & Fabian, A. C. 2005, MNRAS, 358, 211 

Shemmer, O., Uttley, P., Netzer, H., & McHardy, I. M. 2003, MNRAS, 343, 
1341 

Turler, M., Chernyakova, M., Courvoisier, T. J.-L., et al. 2006, A&A, 451, LI 
Tiirler, M., Paltani, S., Courvoisier, T. J.-L., et al. 1999, A&AS, 134, 89 
Turner, M. J. L., Williams, O. R., Courvoisier, T. J. L., et al. 1990, MNRAS, 244, 
310 

Uttley, P., Edelson, R., McHardy, I. M., Peterson, B. M., & Markowitz, A. 2003, 

ApJ, 584, L53 
Walter, R. & Courvoisier, T. J.-L. 1992, A&A, 258, 255 
Winkler, C, Courvoisier, T. J.-L., Di Cocco, G., et al. 2003, A&A, 41 1, LI 
Yaqoob, T. & Serlemitsos, P. 2000, ApJ, 544, L95 

Zdziarski, A. A., Lubiriski, P., Gilfanov, M., & Revnivtsev, M. 2003, MNRAS, 
342, 355 



6. Summary 

In this paper we have presented the results of monitoring of the 
quasar 3C 273 over the period of 2003-2005 over the broad 
energy band, spanning from optical/UV up to soft y-rays with 
XMM-Newton and INTEGRAL. Our monitoring campaign gave 
some surprising results. In particular, in our data we find little 
support for the previously observed correlations between cer- 
tain characteristics of the source, such as the ones found in 
Beppo-SAX data bv lGrandi & Palumbol d2004 ( Fig. [5] ), and in 



GINGA/EXOS AT data bv lWalter & Courvoisier! ( fl992h (Fig. 
Moreover, we note the appearance of a new correlation between 
the UV and the X-ray flux from 3C 273, which was not present 
before. This points to the fact that the source undergoes a secu- 



